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FRACTURE  PROBLEMS  IN  POWER  LAW  VISCOELASTIC  MATERIALS 


It  is  intended  to  continue  the  fracture  work  by  examining  the  Mode 
III  crack  in  an  infinite  sheet.  Moreover,  a  very  general  moduli 
description  will  be  adopted  with  the  Standard  Linear  Solid  (SLS)  and 
the  Power  Law  (PL)  as  special  cases.  The  Riemann-Hilbert  approach 
utilized  in  the  original  work  can  be  as  effectively  utilized  here.  Though 
the  last  step  of  actually  producing  exact,  global  closed-form  expressions 
for  the  field  variables  is  unlikely,  it  is  possible  to  extract  stress 
intensity  factors  and  near-tip  behavior)  It  should  be  noted  that  the 
Weiner-Hopf  approach  is  inadequate  for  this  problem  and  that,  of  course, 
correspondence  principles  are  inapplicable. 

It  is  also  proposed  to  examine  the  debonding  of  an  elastic  plate  on 
a  visoelastic  foundation  as  a  combination  of  lo^  ig,  fracture,  and 
lamination.  The  physical  picture  would  be  ^d  moving  along  a  plate 
and  a  semi-infinite  crack  (whose  tip  might  be  in  front  of  the  load!) 
moving  with  it.  The  viscoelastic  media  will  be  very  general  again. 

Coupled  field  equations  for  the  plate  and  the  half space  must  be  solved. 

The  quantities  of  interest  are  the  position  of  the  crack  tip  (and  its 
stability  in  a  subsequent  analysis),  the  shape  of  the  crack,  and  its 
stress  intensity  factor.  There  are  also  questions  regarding  the 
existence  of  a  solution  for  all  load  speeds. 

Two  papers  resulted  from  the  research  performed  under  this  con¬ 
tract: 

(1)  Nachman,  A.  and  Walton,  J.R. :  Energy  Release  Rate  Calculations 
for  Interface  Edge  Cracks  Based  on  a  Conversation  Integral. 

J.  Intern ' tl.  Solids  and  Structures,  in  press. 

(2)  Nachman,  A.  and  Walton,  J.R.:  Energy  Release  Rate  Cal¬ 
culations  for  an  Interface  Mode  III  Edge  Crack  Based  on  a 


NOTICE:  This  document  was  prepared  under  the  sponsorship  of  the  Air 
Force.  Neither  the  U.S.  Government  nor  any  person  acting  on  behalf  of 
the  U.S.  Government  assumes  any  liability  resulting  from  the  use  of  the 
information  contained  in  this  document. 
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ENERGY  RELEASE  RATS  CALCULATIONS  FOR  INTERFACE  EDGE  CRACKS 
BASED  ON  A  CONSERVATION  INTEGRAL 

A.  Nachman  and  J.  R.  Walton  * 

Introduction 

In  a  recent  paper,  L.  B.  Freund  [1]  demonstrated  the  usefulness  of 
Che  M-integral  conservation  law  in  the  determination  of  stress  intensity 
factors  for  2-diaensional  cracks  in  homogeneous  elastic  bodies  with  certain 
geometric  properties.  For  a  detailed  discussion  of  the  M-integral  and  its 
special  features  chac  permit  the  applications  given  by  Freund,  the  reader 
is  referred  to  that  paper  and  the  references  contained  therein.  Liberal 
use  will  be  made  in  this  paper  of  many  of  the  results  derived  by  Freund. 

What  is  demonstrated  here,  is  that  the  M-integral  can  be  used  to  deter¬ 
mine  the  energy  release  rate  for  certain  interface  cracks  in  much  the  same 
way  as  for  cracks  in  homogeneous  bodies. 

Three  observations  are  needed  for  this  application.  The  first,  due 
to  Smelzer  and  Gurtin  [2],  is  that  the  J-integral  on  a  small  arc  about  an 
interface  crack  is  equal  to  the  energy  release  rate,  as  in  the  homogeneous 
case.  One  difference  for  cracks  in  homogeneous  bodies  is  that  the  J-integral 
may  also  be  related  to  the  stress  intensity  factor,  while  for  interface  cracks, 
Che  J-integral  is  related  to  a  composite  stress  intensity  factor  which 
has  dubious  utility.  (See  Smelzer  and  Gurtin  [2].)  The  second  observation 
involves  the  nature  of  the  far  stress  field  for  bonded  dissimilar  elastic 
wedges,  and  appeals  to  the  analysis  presented  by  Bogy  [3]  for  bonded 
dissimilar  elastic  quarter  planes.  The  third  observation  is  that  the 
integrand  of  the  M-integral  is  continuous  across  bonded  interfaces  lying 
Slong  radial  lines  of  the  chosen  coordinate  system. 

^Department  of  Mathematics ,  Did  Dominion  University,  Norfolk,  VA.  23508. 

*  Department  of  Mathematics,  Texas  ASM  University,  College  Station,  7X  '"943 

Sponsored  by  the  United  States  Army  unfei  Contract  No.  DAAG29-75-C-CC24 . 
Supported  in  part  by  the  United  States  Air  Force  under  AFOSR  Contract 
F4962Q-79-C-0076  and  under  AFOSR  Grant  77-3290. 


4 


Statement  and  Analysis  of  the  Problem 


We  present  here  an  extension  of  the  result  obcained  by  Freund  for  the 
elastic  half-space  with  an  edge  crack  whose  corners  are  sub:ected  to  normal 
and  shear  point  loads.  Specifically,  the  M-incegral  conservation  law  is 
applied  to  the  2-dimensional  plane  strain  problem  of  two  infinite  isotropic 
elastic  wedges  with  opening  angles  w'  and  io",  respectively,  and  with 
different  elastic  properties  (exhibited  through  E’,  u'  and  E",v",  where 
E  and  v  denote  Young's  modulus  and  Poisson's  ratio),  which  are  bonded 
together  along  one  edge  excepc  for  a  crack  of  length  i  extending  from 
the  apex,  and  whose  corners  are  subjected  to  a  system  of  normal  and  shear 
point  loads,  given  by  P',  Q' ,  P"  and  Q".  (See  Fig.  1.)  An  application 
of  the  M-integral  conservation  law  yields  a  relation  among  the  parameters 
P,  Q,  E,  V,  l  and  the  energy  release  rate  of  the  crack.  (Symbols 
wit  out  primes  refer  to  both  materials.)  It  is  clear  that  the  com¬ 
plexity  of  the  problem  precludes  the  determination  of  the  energy  release 
rate  by  first  solving  the  corresponding  boundary  value  problem. 

The  M-integral  is  given  by 


M 


f(WniXi 


TkUk,iXi)ds 


(1) 


where  W  is  the  elastic  energy  density,  ni  is  the  unit  normal  to  C 
(which  we  take  to  be  directed  to  the  right  when  C  is  traversed  in  a  given 
direction)  and  T^  is  the  traction  acting  on  the  material  to  the  left  of 
C.  The  stresses  are  related  to  the  elastic  energy  V  and  the 

strains  •  (u  .  +  «  .)/2  by  2W  ■  a.  .u,  »  where  «  denote  the 

displacements.  The  traction  T^  is  given  by  T^  ■  d^n^.  The  elastic 
body  is  assumed  to  be  in  equilibrium  without  body  forces,  ie  0  -  0, 

»j 

and  che  stress-strain  relation  is  given  by  c. .  m  -r -  .  The  conservation 

ij  ®€ij 


lav  -or  H  is  that  Jf  *  0  whenever  C  is  a  closed  path  surrounding  a 
sin  'onnecced  region  in  che  body. 


i 


Consider  che  integral 


m  -  m'  +  >r 

* 

•  J(WVl  -  Vlc,l*l)ds 


'  ‘•'ft  -  -+  f  <-vvl\  -  rK,A>ds' 


where  the  contours  F'  and  F"  are  as  indicated  in  Fig.  1  and  F  *  F'  J  F" . 

It  follows  easily  from  the  discussion  by  Freund,  that  there  is  zero  contri¬ 
bution  to  M  from  the  parts  of  T  along  the  crack  faces  and  the  outside 

edges  of  the  wedge.  Moreover,  the  value  of  M  on  the  vanishingly  small 

dP 

arc  around  the  crack  tip  is  2— j  ,  where  P(Z)  is  the  potential  energy  of 
the  wedge  with  a  crack  of  length  l.  (See  Rice  [4]  and  Smelzer  and  Curtin  [2].) 
Hence,  on  this  small  arc,  M  is  the  product  of  che  crack  length  and  the 
rate  of  decrease  of  the  energy  with  respect  to  crack  length. 

The  contributions  to  M'  and  M"  from  the  small  arcs  around  the 
corner  points  of  the  wedge  follow  from  a  general  result  derived  by  Freund 
for  an  infinite  elastic  wedge  with  corner  loads.  In  particular,  from 
the  small  corner  arc  in  T'  we  obtain  the  contribution 


-  v2)  f 

e'  ]_y  + 


stnoj1  ii)1  -  simp 


-]  • 


and  from  F" 


-  m"2)  f  Fa" 
E"  J_(i)"  + 


_  ..  <- _  I 

j"  +  u"  -  sintp"  J 


where,  employing  Freund's  notation. 
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F  '  *  -P'cos(u72)  -  Q'sin(u'/2) , 

F  '  *  ?'sin(tj'/2)  -  Q'cos(jj'/2) , 

F  "  *  -P"cos(cu"/2)  -  Q"cos(jj"/2)  , 

Fc"  »  -P"sin(oj"/2)  +  Q"cos(<jj"/2)  . 

On  the  bond  line,  x ^  *  0,  >  l,  the  contributions  to  M  from  F'  and 

F"  cancel,  since,  among  the  stress  and  displacement  components,  only  c.1 
is  discontinuous  across  the  interface,  and  on  that  line  n.x.  =  0  (because 

ii 

the  interface  is  radial)  and  T.u  x.  »  (u  <j  +u  a.  ,)x. . 

K  K,i  1  1,1  11  1,1  Ll  1 

It  remains  to  determine  the  contribution  from  the  large  arc  for  which 
it  suffices  co  know  the  far  field  solution  for  two  wedges  bonded  together 
with  no  crack  and  with  apex  loads  P'  +  P"  and  O'  -  Q".  From  the  analysis 
presented  by  3ogy  [3]  for  two  bonded  quarter  planes,  it  is  clear  that  the 
far  field  is  radial.  More  precisely,  if  the  elastic  fields  are  represented 

with  respect  to  polar  coordinates  (r,0),  then  raQO  and  r a  _  vanish 
for  large  r  uniformly  in  9,  whereas  r^rr  does  not.  Consequently,  on 
the  large  arc  we  may  assume  O^g  and  Ogg  are  zero.  It  is  easy  to  see 
that  in  this  case,  to  satisfy  the  equilibrium  equations  and  the  compara¬ 
bility  equation,  we  must  take  for  the  stress  field 


O'  * 
rr 


•  A-  eos9  -fr  3  sin9 


r9 


99 


0. 


(2) 


The  four  constants  A',  B',  A"  and  B"  may  be  calculated  from  the  four 
conditions 


ug(r,  9  ■*  0-)  -  u^'(r ,9  *  0+) ,  (3) 

uj(r,9  *  0  -)  -  u^(r,9  *  0+) ,  (4) 
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rC  fto" 

ra'  'r,3)eos(*)d9  +  ra"  (r,e)cos(9)d9  -  -(?’  +  ?") , 

,  ri  'o  rr 

-U) 

rO  ru"  i. 

[  ra'  (r,9)sin(-3)d9  +  ra  (r ,9>sin(3)d9  »  (Q’  -  Q”) . 

J  rr  Jn  rr 


C5) 

(6) 


Equations  (3)  and  (4)  assert  continuity  for  the  displacements  on  the  bond 
line,  while  equations  (5)  and  (6)  express  che  equilibrium  of  tractions  on 
the  wedge  — uu '  <  9  <  u",  0  <_  r  <_  R,  where  R  is  the  radius  of  the 
large  arc.  A  simple  calculation  shows  that  (5)  and  (6)  reduce  to 


A' (1-cos  2iu')  -  B'(2u’-sin  2u’)  -  A" (1-cos  2u")  -  B"(2u"-sin  2u")  -  4(Q'  - 
—A '  ( 2ci> '  +  sin  2ui')  +  3' (1-cos  2cu’)  -  A"(2e)"+sin  2cj”)  -  3" (1-cos  2<u")  -  -4(?’ 


Substitution  of  (2)  into  the  polar  fora  of  the  stress-strain  law 
followed  by  Che  application  of  (3)  and  (4)  yields  the  relations 

B’/B"  -  A'  /A"  »  (1  +  a) / (1  -  a),  (9) 

where  a  is  one  of  che  two  Dundurs  bi-material  parameters  given  by 
(See  Bogy  [3].) 


It  is  now  an  easy  matter  to  solve  equations  (7),  (S)  and  (9)  for  A',  3',  A1 
and  S" .  In  particular,  we  obtain 

A'  -  (a(Q'  -  Q")  +  b(?'  +  ?")) (4/d) 

B'  -  (-a(P ’  +  ?")  -  c(Q’  -  Q"))(4/d) 


I")  (7) 

+  ?")•  (3) 
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where  a  *  (1  -  cos  2x’)  -  (y  — (1  -  cos  2u”) 


b  «-<2uj'  -  sin  2cj' )  -  (7-^7)  Ccj"  -  sin  2u>") 

1  +  a 

c  -  -(2u’  +  sin  2oj * )  -  +  sin  2u") 

d  «  a2  -  be. 


A"  and  B"  may  now  be  calculated  from  (9) • 

As  observed  by  Freund,  for  the  stress  state  (2)  ,  the  integrand  of 

M  is 


Wx  n  -  T,u  x  •  -  i  (1  ~  rC2 
11  k  k,i  i  2  E  rr 


Consequently,  the  contribution  to  M  from  the  large  arc  is 


-  \  ^  ~gV’  |  (A'cosS  +  B'sin0)2d6  -  j  (1  ~  ) 


ruj” 


(A"cos0  +  B"sin0)d8 


-<D 


1  (1  -  v')2 

"2  E' 


rO 


(Afcos0  +  B*  sin0)  d0  +  (• 


Lia ,  f 

1  +«  Jn 


(A'cosB  +  B'sin0)d£]  (10) 


-U) 


-  -  -(-  (b(P'  +  P")2  +  c(Q'  -  Q")2  +  2a (P1  +  P")(Q'  -  Q'*))(2/d).  (11) 


Line  (10)  follows  from  (9)  and  Che  observation  chat 


2  2 

^  ~„v — ^  -  (y  —  y)  — — whereas,  line  (11)  is  derived  by  simple  but 
tedious  algebraic  manipulations. 
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Combining  the  contributions  to  M  from  the  large  arc,  the  two  small 


arcs  at  the  apex  and  the  arc  around  the  crack  tip  and  appealing  to  the 
conservation  law,  we  obtain 

2 

z  dP  ,  IL-lJll  I  [_(b(P'  +  p")2  +  c(Q’  -  Q")2  +  2a(P '  +  P"Kq'  -  Q"))(2/d> 

(12) 

+  i(?,2(2u'  -  sin  2m')  +  Q,2(2u’  +  sin  2m')  +  2P'Q'(cos  2m'  -  1)) /(uj,2-sin2y') 

to 

+  i(^-^)(P"2(2a),,-sin  2u")  +  Q"2(2w"  +  sin  2m")  +  2PM0’’(cos  2u"-l)) /(w"2-sin2u")  ] . 
l  l  -  a 


An  important  special  case  of  (12)  is  that  of  two  bonded  dissimilar  quarter 
planes  with  an  edge  interface  crack.  Setting  w'  -  w"  *  r/2  in  (12)  yields 


,dF  (1  -  u'2) 
d!  E' 


[ (P ’  +  PM)2  +  (Q’  -  q")2  -  ^  (p-  +  p") (q1  -  q") ) (l  +  a)/((2a)2~T2) 


+  2((P'2  +  <• 


+  a^p.,2) 


1  -  a- 


+  (Q’2  +  ( 


,l^a 

1  -  a 


)Q") 


-  ^  (P'Qf  + 


A  + 


A  -  a 


a.)P,'Q"))/(m2-4)]. 


(13) 


It  should  be  noted  that  if  P*  -  P"  =  P,  Q’  -  Q"  =  Q  and  a  «  0,  then 
(13)  reduces  to  the  result  obtained  by  Freund  for  identical  quarter  planes. 

Of  course,  when  a  «  0,  the  response  of  the  bonded  quarter  planes  is  the  same 
as  for  a  homogeneous  half-space;  as  described  by  Dundurs  [5],  the  two  quarter 
planes  are  "consonant  in  tension  parallel  to  the  Interface." 

Another  interesting  case  in  (12)  is  when  a  ■  ±1.  Due  to  symmetry  we  consider 
only  a  *  -1  corresponding  to  which  (12)  becomes 
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■  dP  (1 
Ld i  “ 


v'2) 


,.2 


,.2 


jp-  '  [P *  (2w'-sin  :u’)  +  Q'  (2w'+sin  2u')  +  2?'q'C cos  2w,-l))/(w" 


•sin"_  ' 


It  should  noc  be  surprising  that  in  (14)  only  P'  and  Q'  appear,  since 
a  ■  -1  corresponds  to  E"  *  ®. 

It  should  also  be  noted  that  for  P 1  ■  ?"  =  P,  Q'  *  Q"  =  Q  and 
id'  *  at"  *  t/2,  the  level  curves  of  l  in  the  (?,Q)-plane  are  ellipses 
(straight  lines  if  a  ■  0)  centered  at  (0,0).  From  this  we  may  conclude 
as  did  Freund,  that  crack  extension  may  result  when  the  unloading  of  P  and 
Q  occurs  along  certain  paths  in  the  (P,Q) -plane.  Moreover,  this  obviously 
is  the  case  also  in  (14)  with  w'  f  ir/2.  It  is  evident  in  (12), 
that  in  general,  regardless  of  the  values  of  id'  and  w"  and  a,  this 
phenomenon  is  to  be  expected. 

The  authors  wish  to  acknowledge  with  much  appreciation  several  helpful 
and  encouraging  discussions  on  this  problem  with  Professor  L.  3.  Freund. 
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ENERGY  RELEASE  RATE  CALCULATIONS  FOR  AN  INTERFACE 
MODE  III  EDGE  CRACK  BASED  ON  A  CONSERVATION  INTEGRAL 


BY 

A.  Nachman1  and  J.  Walton2 
ABSTRACT 

The  M-integral  is  applied  to  the  calculation  of  energy 
release  rates  for  interface  edge  cracks  of  the  Mode  III  type: 
specifically,  for  an  edge  crack  along  the  interface  between 
two  elastic  wedges  of  different  opening  angles  and  dissimilar 
elastic  properties,  and  that  is  subjected  to  point  loads  at 
the  apex,  a  relation  is  derived  along  the  length  of  the  crack, 
the  energy  release  rate  of  the  crack,  the  applied  loads,  the 
wedge  angles  and  the  material  parameters. 
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In  this  note  we  continue  the  work  of  Cl],  wherein 
the  M-integral  was  used  to  calculate  energy  release 
rates  for  interface  edge  cracks  of  the  Mode  I  and  Mode  II 
types.  However,  here  we  turn  our  attention  to  a  Mode  III 
interface  edge  crack  as  might  arise  in  the  torsion  of  two  welded  rods . 

The  physical  scenario  is  identical  to  Cl]  save  for  the 
loading  at  the  wedge  apexes.  In  this  investigation  we 
impose  point  shear  stresses  (a^)  of  magnitudes  and  P2- 
See  Figure  1.  This  loading  will  generate  only  out-of-plane 
displacements,  u^(x,  y)  .  Moreover,  u^Cx,  '/)  is  a  harmonic 
function  (or,  more  precisely,  (x,  y)  and  u^2^  (x,  y) 

are  each  harmonic  functions  in  their  respective  wedges) . 

We  take  the  same  path  for  the  M- integral  as  in  Cl]. 

As  before,  radial  lines  contribute  nothing,  the  small 

dP 

circle  around  the  crack  tip  contributes  [2],  the  inter¬ 
face  contributions  cancel  and  we  must  only  deduce  the 
contributions  from  the  arcs  near  the  tips  and  from  the 
large  arc. 

Since  the  edges  near  each  tip  are  free,  an  analysis 
of  the  near-tip  field  is  equivalent  to  the  analysis  of  a 
single-material  wedge  subject  to  a  shear-stress  point  load. 
Clearly  the  displacement  field  for  such  a  problem  is 
u<2>  ■  In  (r)  *  3^9  for  the  top  wedce,  say. 
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The  contribution  to  M  from  the  arc  near  the  top 
tip  is  then 


w„'«‘u  9  +  xi  - 1 


where  a 


(2) 


(2)  .  „<21CA<2>c05  8  -  3(2>sin  3]  and  ai?>-  l2> sin  3  * 


31 


32  r 


B  cos  33.  A  simple  calculation  shows  that  the  integral 

equals  u(2)w"  { CB  (2)  ] 2  -  [A(2)]2}.  Since 
2 - 


fW' 


12) 

rarZ  dS  *  -?2cos  ( 


and 


w" 


(2) 


(2)  ^  -P2cos ( j  ) 


ra-  d0  »  -P-sin  (w")  we  have  that  A 
0  32  t  J 


2u(2)^ 


and  B(2)  *  -PjSinCj  )•  Consequently,  the  contribution  to  M 
2\i  ^  2  ^  w" 

from  both  near-tip  arcs  is 


p*  oos  w" 


TU 


p*  oos  w' 


w"  8u  ^  w' 


(1) 


The  contribution  from  the  large  arc  comes  from  the 
analysis,  as  in  Cl],  of  the  far-field  solution  for  two 
wedges  bonded  together  with  no  crack.  Here  the  displacement 
field  in  each  wedge  is  likewise  of  the  form  u^  *  C  ln(r)  +  D9  . 

Since  the  shear  stresses  corresponding  to  the  above 
displacement  field  are 


16 


a31  “  uC^’cos  9  -  -  sin  9) 


a32  *  u  sin  9  +  r  C0S  9) 


(2a, b) 


it  follows  that  continuity  of  u^  and  c^2  at  tiie  interface 
is  assured  when 


CU)  .  c(2) 
u(l)D(l)  =  w(2)d(2) 


Qa,b) 


We  still  need  two  more  equations  and  these  come  from,  as 
in  Cl],  a  force  balance.  Thus 


rd  ■  m  j-w 

ro^;  d9  +  J 

-w'  0 


ra^5  d9  = 


-(Pi  +  P2) 


(4a, b) 


r. 

-w 


ra^  d0  +  JW  ra^5  d8  =  0 


Using  (2a, b)  and  (3a, b)  we  get 


C^cu^  sin  w'  +  u(2)  sin  w"  ] 


+  y (1)  D(1)  Ceos  w  "  -  cos  w?]  =  -  (P^  +  Pj) 


(5a, b) 


C(1)  Cy (1)  (cos  w'  -  1)  +  y(2)(l  -  cos  w")3 


+■  u^D^Csin  w'  +  sin  w"  ]  =  0 


Consequently,  the  contribution  to  M  from  the  large  arc  is 


-U(1)  °  {[C(1)]2  -  CP(1)]2}da  -  y(2)  {[C{2)]2  -  [D(2)]2}d3 

7  -  2  '  0  . 


u(1)w' {cc(1)  ]2  -  CD(1)]2}  -  \  y  (2)w"uc(2)  ]2  -  [d'(2)]2}. 


We  content  ourselves  with  exhibiting  the  particulars 
of  the  case  w',=  w"  =  w.  In  this  case  (=&a,b)-,  together  with 
the  previous  calculations,  serves  to  give 
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We  assume  w/o  l.o.g.  that  u  ^  and  then  some  algebra 

will  show  that  the  level  curves  of  (6)  are  rotated  ellipses. 
Thus,  again,  crack  extension  may  result  when  the  unloading 
of  P^  and  ?2  occurs  along  certain  paths  in  the  (P^,  P2)  plane. 
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